Suitability of hybrid gravitational waveforms for unequal-mass binaries 
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This article studies sufficient accuracy criteria of hybrid post-Newtonian (PN) and numerical 
relativity (NR) waveforms for parameter estimation of strong binary black-hole sources in second- 
generation ground-based gravitational-wave detectors. We investigate equal-mass non-spinning bi- 
naries with a new 33-orbit NR waveform, as well as unequal-mass binaries with mass ratios 2, 3, 4 
and 6. For equal masses, the 33-orbit NR waveform allows us to recover previous results and to 
extend the analysis toward matching at lower frequencies. For unequal masses, the errors between 
different PN approximants increase with mass ratio. Thus, at 3.5PN, hybrids for higher-mass-ratio 
systems would require NR waveforms with many more gravitational-wave (GW) cycles to guarantee 
no adverse impact on parameter estimation. Furthermore, we investigate the potential improve- 
ment in hybrid waveforms that can be expected from 4th order post-Newtonian waveforms, and find 
that knowledge of this 4th post-Newtonian order would significantly improve the accuracy of hybrid 
waveforms. 



I. INTRODUCTION 

Within the next few years, Advanced LIGO (JQ, 
Virgo [|| , and LCGT [f| will likely observe the first grav- 
itational wave signals. The detection rate is expected to 
be somewhere between 0.4 and 1000 detections per year 
for binary black-hole (BBH) systems with masses below 
100 Mq at distances of several gigaparsecs @. Because 
the gravitational wave (GW) signals will be faint com- 
pared to the detector noise, accurate source modeling of 
the predicted GWs will be necessary to detect signals 
and perform parameter estimation on them via matched 
filtering. 

Several methods have been developed to calculate 
the GWs from a BBH system, two of which are post- 
Newtonian (PN) theory and Numerical Relativity (NR). 
The PN expansion Q is a slow-motion, weak-field ap- 
proximation to General Relativity, and it provides an 
accurate description of the inspiral prior to merger, but 
becomes increasingly inaccurate close to merger. NR can 
be used to model the late inspiral, merger, and ringdown 
of a BBH evolution, but at high computational cost 0,0. 

In order to use the best features of each type of model, 
one stitches a long PN inspiral to the front of an NR 
waveform of the late inspiral, merger, and ringdown, thus 
creating a hybrid waveform. The trick is to stitch the 
NR and PN parts together such that the error due to 
higher-order unknown terms in the PN waveform is small, 
while necessitating as few NR orbits as possible to reduce 
the computational cost. This has been explored in our 
previous wo rk llOl for equal-mass, non-spinning binaries, 
and also in (nfllj. 
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Hybrid gravitational waveforms have many uses. They 
play an important role in the creation of phenomenologi- 
cal waveforms [IJ 1111 > which are used in template banks 
for event detection. In addition, hybrids are used within 
the NINJA project fl6l4T8l | in order to test the GW search 
pipeline for ground-based detectors. 

A sufficient criterion to determine whether PN+NR 
hybrid waveforms are suitable for parameter estima- 
tion in the advanced detector era is that a GW model 
ft-modei must be deemed indistinguishable from the true 
waveform h exact within the Advanced LIGO noise spec- 
trum. That is, we place upper bounds on the er- 
ror measure ||<5/i||/||/i|| (to be defined in Sec. Ill Dl) . 
where Sh = h mo dei ~ ^-oxact- This error measure is 
distance-independent as long as cosmological redshifts 
are neglible. 

This paper extends our earlier analysis in several di- 
rections: First, we calculate for equal-mass hy- 
brids constructed from a previously unpublished 33-orbit 
NR waveform in order to confirm results obtained with 
the 15-orbit NR used in prior works. Second, we extend 
the parameter space to hybrid waveforms from unequal- 
mass binaries by calculating this error in the most de- 
manding usage case of parameter estimation. Third, we 
examine the possibility of decreasing hybrid errors if PN 
were known to 4th order. A decrease in PN error would 
mean that more accurate hybrid waveforms could be cre- 
ated with the current length NR waveforms, thus decreas- 
ing the computational cost of generating these templates. 

This paper is organized as follows: Section [TTI reviews 
the methodology of our analysis, including the types of 
PN and NR waveforms, the hybridization procedure and 
how the errors are quantified. Section UHl extends our ear- 
lier results in [lj| and [lj| using a new 33-orbit NR wave- 
form. Section IIVI applies our error analysis to unequal- 
mass binaries. Finally, section IVl estimates the errors one 
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would obtain for PN+NR hybrids if the 4PN terms were 
known. 



II. METHODOLOGY 

A. Post-Newtonian waveforms 

Post-Newtonian (PN) theory presents a slow-motion, 
weak-field approximation to General Relativity in terms 
of expansions of GM/rc 2 and v 2 /c 2 . We use the same PN 
approximants as in our previous work [Toj . Specifically, 
we investigate the properties of TaylorTl, TaylorT2, Tay- 
lorT3, and TaylorT4 (as defined in [ll^Up 3.5PN or- 
der in phase and 3.0PN order in amplitude |22l429| . These 
four approximants differ only in their unknown higher- 
order PN terms. As in [loj . we consider the (2,2) mode 
of the spin- weight s — — 2 spherical-harmonic decompo- 
sition of the gravitational waveform. The amplitude of 
the PN waveform of this mode is always used at 3PN or- 
der, except for Fig. |3] which uses the amplitude to 2.5PN 
order for consistency with earlier work [19( . 
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FIG. 1. The trajectories of the black holes for the 33-orbit 
numerical simulation. The blue line shows the trajectory of 
one black hole and the orange curve shows the trajectory of 
the second black hole. 



C. Hybridization Procedure 



B. Numerical waveforms 

The numerical relativity waveforms were produced 
with the SpEC code [30| , a multi-domain pseudospectral 
code to solve Einstein's equations. We use two simula- 
tions for the equal-mass, non-spinning binaries. One of 
these two simulations covers 15 orbits; it was presented 
in Refs. 0, HH and was already used in the preceding 
study in MacDonald et al. The second equal-mass, 
non-spinning waveform covers 33 inspiral orbits. This 
is a new simulation which is part of a larger, ongoing 
parameter-space study of binary black holes. This wave- 
form was obtained with numerical techniques similar to 
those of [lH . The trajectories of the black holes in the 33- 
orbit simulation are shown in Fig. [1] and the correspond- 
ing gravitational waveform in the top panel of Fig. [2] It is 
more than twice the length of the 15-orbit waveform, and 
will allow us to reduce the GW matching frequency from 
0.038/Af to uj m = 0.025/Af where M = Mi+M 2 with Mi 
and Mi the individual masses of the black holes. In this 
paper, all frequencies u> refer to the gravitational-wave 
frequency that is extracted from the (2, 2) mode of the 
gravitational waveform; at leading order, this frequency 
differs by a factor of two from the orbital frequency often 
used in PN calculations. 

The unequal-mass waveforms of mass ratios 2, 3, 4, 
and 6 were presented in detail in Buchman et al. |32j. 
The simulation with mass ratio 6 is plotted in Fig. [2] it 
covers about 20 orbits. The simulations with mass ratios 
2, 3, and 4 cover are somewhat shorter and cover about 
15 orbits. 



The hybridization procedure used for this investigation 
was the same as in our previous work (see Sec. 3.3 of [Io|): 
The PN waveform, /ipn(^), is time and phase shifted to 
match the NR waveform, /inr(£), and they are smoothly 
joined together in a GW frequency interval centered at 
ojm with width 5ui: 



5uj Suj 

w m - — < LO < CJ m + — . 



(1) 



This translates into a matching interval £ m j n < t < t max 
because the GW frequency continuously increases during 
the inspiral of the binary. As argued in [Io| , we choose 
Slu = 0.1w m because it offers a good compromise of sup- 
pressing residual oscillations in the matching time, while 
still allowing hp^(t) to be matched as closely as possible 
to the beginning of h^m(t). 

The PN waveform depends on a (formal) coalescence 
time, t c , and phase, $ c . These two parameters are de- 
termined by minimizing the GW phase difference in the 
matching interval [i m in,imax] as follows: 



t' c , $' c = argmin / Upv(t; t c , $ c ) - ^W*)) dt, (2) 

where t' c and $>' c are the time and phase parameters for 
the best matching between /ipn(^) and /inr(^), and <f>(t) 
is the phase of the (2,2) mode of the gravitational radi- 
ation. (Because we consider only the (2,2) mode, this 
procedure is identical to time and phase shifting the PN 
waveform until it has best agreement with NR as mea- 
sured by the integral in Eq. [T]) The hybrid waveform is 
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FIG. 2. Some of the waveforms used in this analysis. Shown 
is the real part of the (2,2) mode of the equal-mass, non- 
spinning waveform with 15 orbits [TsL l3ll]. the new 33-orbit 
simulation, and the waveform for a binary with mass ratio 

9 = 6 m. 



then constructed in the form 

hn(t) = T{t)h PN (t; t' c , $' e ) + [1 - Fit)] h NR (t), 
where F(t) is a blending function defined as 
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In this work, we construct all hybrids using the same 
procedure, Eqs. (JT])-(j4]), and we vary only the PN ap- 
proximant and the matching frequency u m . 



D. Quantifying errors 

As in [loj . the error measurement used to deter- 
mine the indistinguishability of two hybrid waveforms 
within a gravitational wave detector's noise spectrum 
was or the noise- weighted inner product of 

the difference Sh between the two hybrids. We denote 
this difference Sh = /ijji — ^H2, where /ihi and /iH2 
are two hybrid waveforms to be compared. The norm 
||<5/i|| 2 = (6h,8h) is defined through the noise-weighted 
inner product 

(9,h)=2j Q df ^ , (5) 

where <?(/) and h(f) are the Fourier transforms of two 
waveforms g(t) and h{t). S n (f) denotes the (one-sided) 
power spectral density, 
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C„(t) , / > 0, (6) 



where C n (r) is the noise correlation matrix for zero- 
mean, stationary noise. As in our previous work, 
we calculate these errors using the Advanced LIGO 
noise curve in its high-power, zero-detuned configuration 
(ZERO _DET_high_P in (H). 

In order to reduce the effects of the Gibbs phenomenon 
in the Fourier transforms, we apply a Planck-taper win- 
dow function (3~H to the time-domain data before com- 
puting the Fourier transform. The error is then mini- 
mized by a time and phase shift of one waveform relative 
to the other. 

Sufficient accuracy of the model waveform is guaran- 
teed if m 
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1/Peff for parameter estimation, 



< 



for event detection. 
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Here, e max is a bound on the fractional signal-to-noise 
ratio (SNR) loss while searching for GW signals. We 
follow the suggestion in [HI and consider e max = 0.005. 
The parameter p e s in Eq. (JT]) represents an effective SNR 
that incorporates a safety factor e < 1 [36j], the impact 
of a network of detectors, and SNR p of the GW event. 
It is defined as 
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N p. 
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where iV is the number of detectors. As in [Io| . we indi- 
cate p e ff = 40 and p e g = 100 in the plots below to cover 
the range of possibilities with a strong GW signal and a 
single to many detectors. In addition, the event-detection 
limit of Eq. (0, V2£max ~ 0.1, can be rewritten in terms 
of Poff = 10. We also indicate this bound in our figures. 



III. RESULTS: EQUAL-MASS, NON-SPINNING 
BINARIES 



As a first consistency check, we perform a PN compar- 
ison similar to Boyle et al. [19j. That work compared 
a 15-orbit NR waveform to the PN approximants Tay- 
lor Tl, T2, T3, and T4. PN and NR waveforms were 
matched at a certain GW frequency Muj m , and then the 
differences in phases between the two waveforms were 
computed as well as the relative amplitude error. We re- 
peat this analysis for the Taylor T3 3.5/2.5PN waveform 
to see if the same behaviour holds for a longer numerical 
waveform. The results are shown in Fig.fSJ As in [l||, we 
matched the NR and PN waveforms together at Muj m = 
0.040, 0.050, and 0.063, but also at the lower frequencies 
0.025 and 0.030. 

The solid lines in Fig. |3] show the phase error obtained 
when using the new 33-orbit NR waveform. The dashed 
lines show the results from [lj|. Agreement between 
these two comparisons is excellent, demonstrating that 
extending the SpEC simulations to larger numbers of or- 
bits yields consistent results with using a shorter wave- 
form. This comparison validates the earlier results [l9j . 
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FIG. 3. Phase error between numerical and post-Newtonian 
waveforms. A TaylorT3 3.5/2.5PN waveform is matched to 
NR at five different GW frequencies u m , and the phase dif- 
ferences between PN and NR are plotted. The solid lines 
show the results obtained using the new 33-orbit numerical 
waveform, with filled diamonds indicating the location of the 
matching point. The thick dashed lines starting 16 orbits 
before merger represent the results achieved in jig ]. 



and our matching procedure, and demonstrates consis- 
tency between the 15- and 33-orbit simulations. 

Continuing the consistency tests between the 15-orbit 
waveform (nil . l3l| and the new 33-orbit waveform pre- 
sented in Sec. Ill Bt Fig. @]shows the error ||<5/i||/||/i|| be- 
tween hybrids created with the two numerical waveforms 
matched at the same hybridization frequency Mu> m = 
0.042 as a function of total mass, as well as the error 
between the high and medium resolutions for each wave- 
form. This figure shows first that the error between hy- 
brids created with the longer and shorter NR waveforms 
is comparable to the numerical error of either simula- 
tion. Numerical errors are much smaller than the error 
bound for p c g = 100, and so we shall therefore disregard 
numerical errors. 

Having established the accuracy of the 33-orbit wave- 
form, let us now assess the quality of PN-NR hybrid 
waveforms at the lower matching frequencies that are 
made accessible by this NR waveform. We choose sev- 
eral matching frequencies in the interval 0.025 < Mu> m < 
0.063, and for each of these, we construct hybrid wave- 
forms with Taylor Tl, T2, T3, and T4. We then com- 
pute pairwise differences ||<5/i||/||/i|| at a fixed total mass 
M = 2OM using the ZERO_DET_HIGH_P LIGO noise- 
curve [33( 1 . The results of this comparison are shown as 
the solid lines in Fig. [5] Differences decrease with de- 
creasing matching frequency (i.e., when switching from 
PN to NR earlier), as one would expect. The thick 
dashed lines represent the same comparison performed 
with the 15-orbit waveform, i.e., the precise data al- 
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FIG. 4. Convergence and consistency tests between the 15- 
and 33- orbit NR waveforms for q=l. The dashed lines show 
truncation error of either waveform, the solid line shows the 
difference between 15- and 33-orbit waveforms. For all three 
curves, hybrids were constructed with the same Taylor ap- 
proximant at the same matching frequency. ||<5/i||/||/i|| was 
then evaluated as a function of mass M using the advanced 
LIGO noise curve. 
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FIG. 5. Error of hybrid waveforms as a function of matching 
frequency for mass ratio q = 1. The solid lines show ||<$/i||/||/i|| 
for hybrids constructed with different PN approximants and 
with the new 33-orbit NR waveform. The dashed lines show 
the previous results which use a 15-orbit NR waveform. These 
errors are all for binary black-hole systems with a total mass 
of 20M o . 



ready presented in Ref. |i0| . The comparisons agree for 
Mu m > 0.038, where both NR waveforms are available. 
Furthermore, the long waveform continues the trend set 
by the earlier, shorter waveform without surprises. This 
indicates that extrapolating errors in our analyses to 
lower matching frequencies is indeed valid. It also in- 
dicates that there is no unexpected behaviour of PN 
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FIG. 6. The relationship between the number of orbits before 
merger and the gravitational wave frequency for mass ratios 
9=1 and q = 6. 



or NR in the newly accessible GW frequency interval 
0.025 < Mu < 0.038, with PN and NR converging to- 
ward each other. 



IV. UNEQUAL-MASS BLACK HOLES 

Now let us consider binary black-hole systems with 
mass ratio q = Mi/M 2 = 2, 3, 4, 6. As q increases, the in- 
spiral proceeds more slowly, in proportion to the symmet- 
ric mass ratio r\ = M1M2/M 2 (see e.g., post-Newtonian 
expansions Q). This is illustrated by Fig. [6l which shows 
the number of inspiral orbits to merger as a function of 
gravitational wave frequency This figure was created 
by plotting the number of orbits prior to the maximum 
amplitude of hybrid waveforms matched at the earliest 
possible GW frequency (uj m — 0.025 for q — 1 and 0.042 
for 9 = 6) with a Taylor T4 waveform against uj. 

Starting from the same GW frequency (e.g., Mui = 
0.046), the q = 6 binary proceeds through roughly twice 
as many orbits. Conversely, the same number of orbits 
to merger (e.g., N = 12) occurs at higher frequency for 
<7 = 6 than for q = 1. The dashed lines indicate a match- 
ing frequency of Mu> m = 0.046 and 12 orbits; we will 
use these two reference values in subsequent comparisons. 
Because unequal-mass binaries spend more orbits in the 
strong field regime, we would expect that with increasing 
mass ratio q, hybrids need to be matched a larger number 
of orbits before merger than for q = 1. The next sections 
will quantify this expectation. 



lated phase difference between PN and NR waveforms. 
Following Hannam et al. [37||, we match PN and NR at 
Mu> m = 0.1 and then calculate their phase difference at 
a certain time before this matching point. 

We perform this computation for TaylorTl, TaylorT2, 
TaylorT3, and TaylorT4. (The TaylorT3 waveform was 
only calculated to 3.0PN order in phase since the 3.5PN 
order waveform does not reach Mui = 0.1.) The results 
are shown in Fig. [7J The top panel shows the phase 
difference between PN and NR at 8 GW cycles before 
Mu) m = 0.1 (equivalent to Fig. 8 of (37j). and the bot- 
tom panel shows the phase difference at Muj = 0.05 for 
PN and NR waveforms matched in phase at Muj m = 0.1. 
Both panels show similar trends. In the lower panel, the 
phase difference increases more rapidly with increasing 
mass ratio because the number of orbits within the com- 
parison increases with increasing mass ratio (7.4 orbits 
for q = 1 versus 12.2 orbits for q = 6). 

Our findings disagree with the results of Hannam 
et al. [37| . despite following the identical comparison 
protocol. Figure 8 of [13] shows a roughly constant 
phase difference between TaylorTl and numerical sim- 
ulation of ~ 0.5 radians, for all considered mass ratios 
(q = 1, 2, 3, 4), whereas we find a steadily declining phase 
difference reaching zero near q ~ 5. Similar trends hold 
for TaylorT4: Ref. [13] reports phase differences ~ —0.1 
radians for mass ratios q = 1, 2, 3,4. We see a very small 
phase difference at mass ratio q = 1, and steadily in- 
creasing phase differences at larger mass ratios. We be- 
lieve that these discrepancies are cause d by the older, less 
accurate NR simulations used in Ref. [371 ] . 

Figure [7] shows that one cannot assume that any sin- 
gle PN approximant will be suitable for all of parameter 
space. Phase errors between PN and NR change dra- 
matically with increasing mass ratio. The agreement be- 
tween TaylorT4 and NR waveforms in the case of equal- 
mass non-spinning binaries, for example, is purely co- 
incidental. As the mass ratio of the binary increases, 
this phase error becomes much larger for Taylor T4, and 
in fact becomes smaller for TaylorTl and TaylorT2. Be- 
cause the Taylor approximants differ only in higher order 
post-Newtonian terms, the spread between them can be 
taken as a measure of the post-Newtonian truncation er- 
ror. For the comparison in the top panel, this indicates 
a post-Newtonian truncation error of ~ 0.5 radians for 
5=1 increasing to ~ 1 radian for q = 6. Within this 
(admittedly large) truncation error, all four Taylor ap- 
proximants are consistent with the numerical data (i.e., 
consistent with zero phase difference). 



B. Hybrid errors 



A. Phase errors between PN and NR 

Let us start with investigating the phase difference be- 
tween NR and PN with measures that are independent 
of the LIGO noise spectrum. We calculate the accumu- 



We now repeat the analysis of Fig. [5] for binaries with 
higher mass ratios: We hybridize TaylorT[l,2,3,4] at sev- 
eral matching frequencies uj m . At each u m , we compute 
differences ||<5/i|j/||/i|| between all six pairs of PN approx- 
imants and plot these differences as a function of to m . 
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FIG. 7. Accumulated phase difference between numerical and 
post-Newtonian waveforms for Taylor Tl, T2, T3, and T4. 
PN and NR are matched at GW frequency Mu: m = 0.1, phase 
differences are then computed eight GW cycles earlier (top 
panel) or at GW frequency Mui — 0.05 (bottom panel). The 
Taylor T3 waveform used in this comparison is 3.0 PN order 
in phase. 



The results arc shown in Fig. [5] for mass ratio q = 6 
and for total masses 1OM and 40M Q . M = 1OM 
represents a binary with component masses 1.67M and 
8.33A/ , which we shall take as an approximation of a 
black-hole-neutron-star (BH-NS) binary. M = 4OM 
results in component masses 6.67M and 33.337M , a 
BBH system where the smaller black hole is consistent 
with known black-hole masses. We chose these masses 
because they are more astrophysically probable for this 
mass ratio than the total mass of 2OM that we used for 
q = 1. The coloured lines represent the differences be- 
tween hybrids as just described. To ease the comparison 
with the q = 1 results of Fig. \5[ we duplicate those data 
into Fig. [8] as the grey lines in the background. 

One notices immediately two differences between q = 1 
and q = 6: (1) At the same matching frequency, q = 6 
results in larger differences. This might be caused by the 
larger number of orbits that the q = 6 binary spends 
at high frequency cf. Fig. [Tj] (2) The q = 6 compari- 
son covers only comparatively high matching frequencies 
Mu> m > 0.042, whereas the q = 1 comparison reaches 
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FIG. 8. Error of hybrid waveforms as a function of match- 
ing frequency for mass ratio q = 6. Plotted is the value of 
||<5/i||/||/i|| at IOMq and 40Mq as a function of u: m for a bi- 
nary black-hole system with a mass ratio q = 6. The data 
from Fig. [5] are plotted in grey in the background for refer- 
ence. 



much lower frequencies. This originates in the slower in- 
spiral of higher mass ratios (i.e., longer time to merger 
from the same starting frequency) and the higher com- 
putational cost of high mass ratio simulations. 

The difference between the TaylorTl and TaylorT2 hy- 
brids is particularly small in Fig. [8] This is consistent 
with Fig. [7j where for q = 6, the differences between PN 
and NR arc similar for these two approximants. 

Let us now investigate the dependence on q in more 
detail. We compute the differences ||<5/i||/||/i|| for all mass 
ratios at the same matching frequency Moj m = 0.046 
(the data for q = 1 and 6 can be read off of Figs. [5] 
and HI respectively) . The upper panel of Fig. [9] shows 
these differences as a function of q. Once again, it is 
evident that differences increase with mass ratio. 

Length requirements for NR waveforms are often 
phrased in the convenient unit "number of NR orbits" . 
To place this notion in context, we match PN+NR 12 
orbits before merger, compute differences ||<5/i||/||/i|| , and 
plot these in the lower panel of Fig. [SJ When match- 
ing a fixed number of orbits before merger, ||<5/i||/||/i|| in- 
creases even more steeply with q, because the matching 
frequency increases with q, cf. Fig. [5] 

Both Figs. [7] and IH1 show that the error between hybrids 
using different PN approximants increases with mass ra- 
tio. Therefore, with increasing mass ratio binaries, NR 
waveforms will have to be longer to create hybrids of sim- 
ilar quality. For mass ratio 6, one might estimate from 
Fig. [8] that total mass M = 4OM requires a matching 
frequency of Mu> m w 0.03. The matching interval Eq. JTJ) 
would then extend to a lower frequency Muj = 0.0285 
requiring about 40 orbits covered by NR. The needed 
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FIG. 9. Error of hybrid waveforms as a function of mass ratio 
q. PN and NR is matched at frequency Mui m = 0.046 (top 
panel) or 12 orbits before merger (bottom panel). Pairwise 
differences between hybrids in either category are computed 
at total mass M = 20M o . 



matching frequency (and thus the number of orbits) de- 
pends on the total mass considered; for M = 10A/ Q , 
convergence of the errors with decreasing Muj m is not 
yet apparent (see the 10M Q curves in Figure [8]), so the 
NR waveform is too short to even estimate how long it 
should be. This indicates that BH-NS systems may very 
well place the most stringent requirements on NR simula- 
tions. (A proper treatment of BH-NS systems, of course, 
will also require to simulate the neutron star directly, 
including its tides and other effects arising from micro- 
physics. Such a simulation would be yet more challenging 
than our approach of using the easier BBH system as a 
proxy.) 

It is no surprise that the error for hybrids matched 
at a certain Muj m or at a fixed number of orbits before 
merger would increase with mass ratio since the number 
of orbits spent in the strong field regime increases with 
mass ratio. 
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FIG. 10. The magnitude of the PN coefficients for Taylor T4 
as a function of PN order. Shown are different mass ratios 
and the test-mass limit. The straight black line shows our 
assumed coefficients for estimating the accuracy of 4PN and 
4.5PN hybrids. 



V. HIGHER-ORDER POST-NEWTONIAN 

The primary source of error lies in the truncation error 
of the PN approximants. Work is currently being done 
to calculate PN to 4th order [H, HU , therefore, an inter- 
esting question arises: how much might higher PN orders 
improve the accuracy of PN+NR hybrid waveforms? To 
address this question we consider the TaylorT4 approxi- 
mant, where the phase evolution is determined by a single 
Taylor series (ioj : 



dx 
~dt 



64c 3 ^ 5 I /743 11 \ t 3/2 
5GM ~~ ( 336 + T x + A-KX 3/2 
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A,x 5 



(9) 



Here, v — M1M2/M 2 denotes the symmetric mass ratio, 
7 is Euler's constant, c is the speed of light, G is the 
gravitational constant, and x = v 2 /c 2 . 

In Eq. [5J we have included the terms A 4 , A4.5, and 
A5 at 4th through 5th PN order. These terms are cur- 
rently unknown, but trends in the known coefficients can 
be used to approximate the magnitude of the unknown 
coefficients. Figure [TU] plots the known coefficients Ak 
vs. PN order k for mass ratios q = 1, 2, 3, 4, 6. The figure 
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also shows the coefficients for the test-mass limit v 0, 
which are known up to 5.5-th PN order [4l|. The co- 
efficients At behave rather erratically, but overall, they 
seem to be exponentially increasing, as indicated by the 
black line with filled circles. This black line represents 
A4 = 300, with a doubling of the coefficients with each 
increase in PN order (and halving with each decrease). 
Therefore, we shall take the unknown coefficients to be 
A± = 300, A4.5 = 600 and A 5 = 1200. These values 
are of course not the correct ones, but are indicative of 
the expected magnitudes of these coefficients. Hence, we 
will not be able to compute the correct 4PN (and higher) 
PN waveforms, but merely estimate the errors in those 
waveforms, were the coefficients known. 

Assuming A4 = 300, we can estimate the truncation er- 
ror of 3.5PN TaylorT4 by computing ||<fa||/||ft|| between 
the standard TaylorT4 (with A k = for k > 4) and 
a modified TaylorT4 with A4 = 300. This comparison 
is shown as the blue curve in Fig. [TT] which is labeled 
"3.5PN" (the label indicates that this is an estimate of 
the error of 3.5PN order). This new estimate of the 
3.5PN truncation error should of course be consistent 
with our earlier estimates shown in Fig. [5] To demon- 
strate this consistency, we include the data of Fig. [5] as 
the greyed out lines in the background of Fig. [IT] Indeed, 
the new estimate ("3.5PN") follows closely the trends of 
the more exhaustive study, lending confidence in this ap- 
proach. 

Repeating this procedure at next higher PN order will 
now result in an error estimate of 4PN (were it known). 
Thus we compare TaylorT4 hybrids with and without 
the A4.5 = 600 term. This results in the red line labelled 
"4PN" in Fig. [TTJ We can go one PN order further, 
and include an A§ term (resulting in the line labelled 
"4.5PN"). We can also remove the (known) 3.5 PN term, 
to estimate the 3PN truncation error if 3.5PN were not 
known (the curve labelled "3PN"; this curve compares 
A 3 . 5 = with A 3 . 5 = 150). 

We repeated this analysis for the mass ratios q = 
2,3,4, and 6. The results are similar to the q — 1 cal- 
culation of Fig. [TTJ as an example, Fig. [12] presents the 
analogous calculation for q = 6. In Fig. [TT] and [12] 
a clear pattern emerges: Each additional PN order re- 
duces |j<5/i||/||/i|| by approximately a factor of ~ 2. When 
matched at low frequency Mu> m (where PN is more accu- 
rate) the reduction in error is somewhat faster than when 
matching at high frequency. To make these statements 
quantitative, Figure [T3] plots the ratio of the "3.5PN" 
and the "4PN" curves in Fig. [TT] and [12] It also shows 
data for the remaining mass ratios (q = 2, 3, 4). This ra- 
tio is 0.5 at Muj m rj 0.04 and drops to 0.42 at the lowest 
accessible matching frequency Mu m = 0.025. The gain 
of higher-order PN is approximately independent of mass 
ratio, except for high matching frequencies and high mass 
ratios; in this regime errors are so large that the asymp- 
totic trends for small errors/early matching arc masked. 
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FIG. 11. Estimated error of hybrid waveforms if PN up to the 
specified order were known. Shown is the equal-mass case. 
The line "3.5PN" represents the error in the presently avail- 
able PN waveforms versus an estimated 4PN term; see the 
text for details. The grey lines show the more exhaustive 
analysis from Fig. [S] which is consistent with the "3.5PN" 
line obtained with our alternative estimation procedure. 
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FIG. 12. Estimated error of hybrid waveforms if PN up to the 
specified order were known. Shown is the mass ratio q = 6 
case. The line "3.5PN" represents the error in the presently 
available PN waveforms. The grey lines show the more ex- 
haustive analysis from Fig. [8] but at a total mass of 20M© • 
Compare to Fig. [TTJ 



A. 3PN Hybrids 

To further quantify the reliability of our estimate of 
4PN improvements, let us apply our analysis to 3PN 
waveforms. The 3PN truncation error estimate based on 
modifying the 3.5PN A3. 5 coefficient is already plotted in 
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FIG. 13. Estimated reduction in error of hybrid wave- 
forms constructed using 4PN relative to those constructed 
with 3.5PN. This reduction is shown as a function of match- 
ing frequency for different mass ratios. The q = 1 line is the 
ratio between the "4PN" and "3.5PN" curves in Fig. [11] 
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FIG. 14. Error analysis of PN-NR hybrids using only 3PN in- 
formation. The solid and dashed blue lines repeat the analysis 
of Fig. [5] at 3PN order. The dotted red lines show ||<5ft||/||ft|| 
for hybrids using 3PN waveforms compared to one using Tay- 
lor T4 at 3.5PN (the circle is Tl, the star T2, the square T3 
and the downward triangle is T4). For reference, the data of 
Fig. [5] are included in grey in the background. 



Figs. [Til and IT2"1 What remains is the equivalent of Fig. [5] 
at 3PN order: We prepare TaylorT[l,2,3,4] hybrids at 
3PN order, and compute their pairwise differences. This 
results in the thick solid blue lines of Fig. [14] Compared 
to the 3.5PN comparison of Fig. [5] (replicated in the grey 
lines in Fig. fT4"]) . the error of 3PN is indeed somewhat 
larger than for 3.5PN. The change in error 3PN vs. 3.5PN 
is smaller than the factor of 2 we would have predicted 
from Fig. [TTJ presumably because the actual A 3 5 coeffi- 
cient is smaller than the assumed A3. 5 = 150 of Fig. QJJ 
Fig. [JJ] offers another cautionary lesson: Consider 
the differences between hybrids using the 3PN Taylor 
T2,T3,T4 approximants. There are three such differ- 



ences, plotted with blue dashed lines in Fig. [JJ] The 
differences between these three hybrids are surprisingly 
small, with ||<5/i||/|j/i| < 0.07, even when hybridized at 
the very large matching frequency Mtu m = 0.065. If one 
had used only these three approximants, and if only 3PN 
were available, one would have concluded that the result- 
ing hybrids are good for Advanced LIGO purposes, and 
that short NR simulations are sufficient. 

But this conclusion would have been entirely incor- 
rect! The dotted red lines in Fig. [Ml show the difference 
between 3PN hybrids with the 3.5PN TaylorT4 hybrid. 
These differences are an order of magnitude larger than 
the 3PN internal TaylorT2, T3, and T4 differences. At 
3PN, TaylorT2, T3, and T4 are very similar, and they are 
all three off by the same large amount. The 3PN Tay- 
lorTl comparisons reveal this effect: As shown by the 
solid blue lines in Fig. [JJ] differences between TaylorTl 
hybrids and the other three hybrids are all very large. 
As this analysis demonstrates, these large differences are 
not caused by a "deficiency" of TaylorTl, but rather by 
coincidentally similar deviations from the true waveform 
in the TaylorT2,3,4 approximants. The same effect at 
3.5PN order is apparent in Fig. 5: At low matching fre- 
quencies, the pairwise differences between TaylorT2,3,4 
are significantly smaller than the differences relative to 
TaylorTl. Therefore, it is important to investigate many 
different approximants, and not to discount lone outliers 
as in Fig. [JJ] 



VI. DISCUSSION 

In this paper, we have presented an analysis of errors 
that affect hybrid gravitational waveforms for a range 
of mass ratios to assess their suitability for parameter 
estimation with Advanced LIGO. We have also estimated 
by how much these errors would be reduced if PN were 
known to a higher order. 

In the case of equal-mass, non-spinning binaries, we 
have found that the results obtained with the most re- 
cent, 33-orbit numerical waveform are consistent with 
previous results obtained using the older, shorter 15-orbit 
simulation. The errors between hybrids using these two 
NR waveforms are small enough for parameter estimation 
for sources with p e g < 100. In addition, when compared 
to PN, these two NR waveforms yield similar results. The 
results in Fig. [5] show that the results using the 15-orbit 
waveform can be reproduced and extended. 

Expanding the parameter space to unequal-mass bina- 
ries, we have found that PN errors grow with increased 
mass ratio as observed by pTl4l4| . Phase differences be- 
tween PN and NR vary strongly with mass ratio and PN 
approximant. Therefore, no single PN waveform is ap- 
propriate for all of parameter space. For example, despite 
the fact that Taylor T4 matches remarkably well with NR 
in the case of equal-mass, non-spinning binaries, this is 
no longer true as the mass ratio between the black holes 
in the binary system increases. It is only a coincidence 
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that TaylorT4 and NR are so similar in this very unique 
configuration. In fact, for q = 6, TaylorTl and TaylorT2 
agree very well with NR, whereas TaylorT4 does not, cf. 

Fig.m 

When evaluating |j<5/ij|/||/ij| for hybrids which use dif- 
ferent PN approximants, it becomes clear that higher 
mass ratio binaries will require increasingly more NR or- 
bits to reach similar accuracy. This becomes problematic 
for high mass ratio binaries, since they are more compu- 
tationally expensive because it requires more numerical 
resolution to resolve the smaller black hole. We also note 
the recent analysis [42j of an equal-mass aligned spin 
BBH simulation with very large spins. Ref. |42| found 
that the simulation covering ~ 25 orbits was of insuffi- 
cient length to reliably hybridize with PN. 

Knowledge of higher-order PN waveforms will signifi- 
cantly improve the quality of PN+NR hybrid waveforms. 
As demonstrated in Figs. [TT] and Q21 any further ad- 
ditional PN order should decrease errors in our mea- 
sure ||<5/i||/||/i|| by roughly a factor of ~ 2, for the same 
length of the NR waveform. Because mismatches are 
proportional to the square of ||(5/i||/||/i|| this will reduce 
mismatches by a factor ~ 4. For a specified accuracy, 
the knowledge of a higher-order PN expansion would al- 
low one to shorten the length of NR simulations. For 
PcH = 10, Fig. [TT] indicates that the matching frequency 
could be raised from Mu> m = 0.027 to Muj m = 0.038, 
thus approximately reducing the temporal length of the 
NR simulation by a factor of ~ 2. This would be a sub- 
stantial saving for future NR simulations. We emphasize 
that these estimates depend on our assumption of the 
approximate magnitude of the unknown PN coefficients, 
as discussed in the context of Fig. [TU] 

This analysis has also provided us with an important 
cautionary tale: one cannot ignore outliers when compar- 
ing many different PN approximants to each other. This 
is illustrated in the 3PN case in Fig. [TT] where the error 
between Taylor T2, T3, and T4 hybrids is very small, 
but when compared to Taylor Tl hybrids, the error is 
much higher. Thus, it is very important, when doing 
this type of analysis, to consider as many different PN 
approximants as possible. 

This work could be extended in a few ways. First of all, 
it would be interesting to further extend the parameter 
space to spinning binaries, and also to the most general 
case of processing binaries. Some work has been done to 



create hybrids for precessing binaries in [43l - l46l | . but er- 
ror estimates on this type of hybrid are still a ways off. It 
would also be useful to extend the error analysis of GW 
modes to modes different from just the (2,2) mode, be- 
cause these other modes become increasingly important 
with higher mass ratio (e.g.. |3 2fl) and with precession of 
the orbital plane (e.g., |40l. |471|). 

It would also be useful to refine our error limit in Eq. [7] 
to incorporate the effects of a network of detectors in a 
more effective way. It might be worthwhile to have an 
error limit which depends on the total mass of the binary 
as in [l3j , or to consider some better way of finding an 
upper bound on ||<$/i||/||/i||. An important consideration 
is the fact that our error criterion is sufficient but not nec- 
essary to the suitability of these hybrid waveforms. For 
instance, investigations into detection efficiency of hybrid 
waveforms flil |4a | indicate that hybrid waveforms can be 
perfectly usable, even when failing the indistinguishabil- 
ity test by a wide margin. 
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